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ABSTRACT 

Using the Hubble Space Telescope ACS imaging of the GOODS North and South helds during Cycles 11, 
12, and 13, we derive empirical constraints on the delay-time distribution function for type la supernovae. 
We extend our previous analysis to the three-year sample of 56 SNe la over the range 0.2 < z < 1.8, 
using a Markov chain Monte Carlo to determine the best-fit unimodal delay-time distribution function. 
The test, which ultimately compares the star formation rate density history to the unbinned volumetric 
SN la rate history from the GOODS/HST-SN survey, reveals a SN la delay-time distribution that is 
tightly confined to 3 — 4 Gyrs (to > 95% confidence). This result is difficult to resolve with any intrinsic 
delay-time distribution function (bimodal or otherwise), in which a substantial fraction (e.g., > 10%) of 
events are "prompt" , requiring less than approximately 1 Gyr to develop from formation to explosion. 
The result is, however, strongly motivated by the decline in the number of SNe la at z > 1.2. Sub-samples 
of the HST-SN data confined to lower redshifts (z < 1) show plausible delay-time distributions that are 
dominated by prompt events, which is more consistent with results from low-redshift supernova samples 
and supernova host galaxy properties. Scenarios in which a substantial fraction of z > 1.2 supernovae are 
extraordinarily obscured by dust may partly explain the differences in low-z and high-z results. Other 
possible resolutions may include environmental dependencies (such as gas-phase metallicity) that affect 
the progenitor mechanism efficiency, especially in the early universe. 

Subject headings: supernovae: general- Accepted to the Astrophysical Journal 



1. introduction 

The progenitor systems and mechanisms responsible for 
type la supernovae (SNe la) remains unresolved. The 
general consensus is that SNe la stem from C+O white 
dwarf stars (WD) that accrete mass until they exceed 
the the electron degeneracy pressure limit in their cores, 
marked by the Chandrasekhar mass limit of approximately 
1.4 Mq. But the details of how the additional mass is ac- 
cumulated, or specifically what the donor mass source is, 
remains largely ambiguous. Very broadly, progenitor mod- 
els are categorized as either involving mass accretion from 
companion stars (typically red-giant stars) in single degen- 
erate (SD) scenarios, or pairs of WDs that merge through 
coalescence or collisions in double d egenerate (DD) sce- 
narios (see lYungelson fc Liviol |2000| for a review). The 
various detailed modeling of these scenarios have thus far 
provided adequate agreement with the observed spectra 
and light-curves of SN la events, largely due to gross re- 
quir ement radioactive Fe- peak material to power th e event 
(see iRosswog et al.l 120091 and iRaskin et al.l l2009bl for re- 
cent examples involving collisional WD mergers). There 
remains much uncertainty as to which scenarios are actu- 
ally employed to make SNe la. 

Unlike the progenitors of core-collapse supernovae that 
are now directly found through deep archival imaging at 
a rate of a few per year, SNe la are much rarer (by about 
a factor of ten in typical low-z galaxies) and their pro- 
genitors are much fainter (by a factor of several million), 



making it extraordinarily unlikely that these progenitors 
will be similarly resolved in the near future. Nonetheless, 
meaningful constraints on SN la progenitors can be drawn 
from the resolvable hosts environments of these events. 
Age limits on the stellar population, the rate of formation 
of new stars, and the range of chemical enrichment in the 
environment of the event all provide implied but impor- 
tant constraints on the nature of progenitor systems — the 
types of stars involved, and the physical mechanisms em- 
ployed to result in these luminous explosions. It is there- 
fore expected that correlations drawn from these environ- 
mental characteristics and characteristics of SNe la (e.g., 
event luminosity, or event production) will eventually illu- 
minate how SNe la are formed. 

However presently, the analysis of SN la rates in high 
and low redshift galaxies show inconsistent results on the 
implied progenitor mechanisms responsible for producing 
these important cosmological tools. At its heart, the dis- 
crepancy hinges on two important factors: (1) the incuba- 
tion time of SNe la (commonly called the "delay time"), 
or the time required for a progenitor system to develop 
into an explosion from a single episode of star-formation; 
and (2) the metallicity of the progenitor at formation, 
and its impact on either production efficiency or event 
luminosity. While these are conceptually measurable fac- 
tors i n low-redshift (z < . 1) galaxies, attem pts to do so 
(e.g., iGallagher etall 120081: iHowell et al]|2009ft have been 
muddled by two degenerate effects: (1) population age, 
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which steadily increases the range of metallicity within a 
given environment, and (2) rate of active star formation, 
which mix-up the time between events and progenitor for- 
mation. 

It is expected that at high redshifts (z > 1) these con- 
fusion effects should become less severe as the age of the 
parent population is limited by the age of the universe at 
z > 1 to be less than 6 Gyr old. Moreover, any substantial 
bulk delay-time would limit the SN la rate in the highest 
redshift regimes as the Universe would not be old enough 
to produce them. Thus observations in the highest red- 
shifts could provide the best leverage in discerning the ages 
of SN la progenitors. However the implied delay-time dis- 
tribution from an investigation of the first year (Cycle 11) 
imaging of the GOODS North and South fields with the 
Hubb le Space Telescope (HST) and ACS (|Strolger et all 
2004) were on average very long (3 to 4 Gyr), and in- 
consistent with the relatively short times predicted from 
binar y star evolutionary models (cf. IHan fc Podsiadlowskil 
2008j). They were also difficult to reconcile with several 
observations of SNe la in low-z galaxies (cf. lHamuv et al.l 
2000 tlMannucci et al ll2005t IScannapieco fc Bildstenll2005t 
Gallagher et al. 2008) that suggest there are at least two 



mechanisms for SN la production, one which could re- 
quire a few Gyr incubation, and a second which is much 
prompter, requiring only a few 100 Myr of incubation. Ad- 
ditionally there was seemingly little support for the im- 
plied trend of an increasingly prompt-dominated mecha- 
nism as obs ervations are pushe d to higher redshifts from 
the ground (jHowell et al.ll2007D . 

With just 25 events from the Cycle 11 data, and few 
events at z > 1.2, it is possible the sample provided too 
little statistical certainty to ascertain the inherent delay 
time distribution function. In addition, the tested models 
were simplistic and did little to address more than the av- 
erage incubation times for SN la progenitors. More flexible 
functionality to the model tests would do better to assess 
the delay time distribution. 

The addition of the Cycle 12 & 13 HST-SN data pro- 
vides an opportunity to revisit the best delay-time model 
analysis on a larger statistical sample with a more com- 
plete analysis. Here we present th e extension of the in- 
vestigation of IStrolger et al.l ([2004), applied to the com- 
plete three-year (Cycles 11, 12, and 13) sample of 56 
SNe la. We compare the observed redshifts of the SNe la 
(in the range 0.2 < z < 1.8) to event redshift distribu- 
tions forecasted from model delay-time distribution func- 
ti ons, assuming the star formation rate density model used 
in IStrolger et all (|2004l ). Through a Markov chain Monte 
Carlo test, we determine the most likely model delay-time 
function, empirically implying the bulk distribution of in- 
cubation times of SNe la. In Sj5] we discuss measurements 
of the SN la rate history, and its relation to the rate of 
star formation and the delay times of SNe la. In S}3] we 
describe the model delay-time tests. In Sj3] we show the 
highest likelihood delay-time distribution model from the 
HST-SN data. And in 2]we discuss the results in compari- 
son to other determinations of the ages of SN la progenitor 
systems. 



2. THE SN IA RATE HISTORY 

The various SD and DD scenarios are expected to 
involve very diff erent characteristic development time s 
and distributions fYungelson fc Livioll2000l:lGreggioll2005l ). 
The rate of SNe la at any given epoch should therefore 
mimic the stellar birthrate at earlier cosmic epochs, shifted 
and convolved with the delay-time distribution function, 
$(r), that is set by the progenitor mechanism scenario. 
Here r is the elapsed time between the progenitor system 
formation and explosion as a SN la event. The volumetric 
SN la rate (Ri a , in units of events yr™ 1 Mpc~ 3 h 3 ) at any 
specific cosmic age can therefore be described by, 



Rla(i) = / $(*-*')/>*(*') 

Jtg 



Ai a / £(M)dM 



M 



dt\ (1) 



where the stellar birth rate is the combination of the 
star formation rate density [hereafter p+(t), in Mq yr -1 
Mpc -3 h], and the initial mass function [hereafter £(M), 
in Mq 1 ]. This equation is easily translated to redshift 
space by defining t as the age of the universe at redshift 
z, and setting to as the age of the universe when the first 
stars were born, corresponding to z w 10. It follows that 
J p*(t) dt = J p±(z) dz. The $(t) distribution function is 
normalized so that the integral sum of the function over 
all incubation times is unity. 

Although there i s evidence that suggests £ (M) m ay vary 
with cosmic time (iDavel 120081: | van Dokkumll2008l) . these 
variations are notably small (with small changes in power- 
law slope) in the < 8 Mq mass range for SN la progen- 
itor system stars. For the purposes of this investigation, 
£(M) is assumed to be more or less invariant with cos- 
mic time. It is also assumed that the SN la mechanism, 
although not ubiquitously affecting all stars in the pro- 
genitor mass range, affects essentially the same fraction of 
stars in £(M) in all cosmic epochs, thus allowing Ai a (t) = 
constant. With these preliminary assumptions, we define: 



Ai a / £(Af ) dM 



M 



(2) 



where e is a constant that describes the number of SNe la 
produced per Mq formed, or an efficiency for the stellar 
progenitor population in actually producing SN la events, 
as presumably not all WDs become SNe la. The rate of 
SNe la can therefore be simplified as: 



Ria(i) =e f $(t-t') Pit (t')dt', 



(3) 



However, recent results from the Sloan Digital Sky Survey- II IDildav 
underway and in the near future. 



where the shape of the rate distribution with time is only 
dependent on star formation rate density history and the 
distribution of delay times. 

With the advancing knowledge of p*(z) and Ri a (z), it is 
plausible to attempt to "deconvolve" Equation [3] to learn 
about $(r), and thus constrain the nature of SN la pro- 
genitor systems, which is the goal of this investigation. 
To this end, we are fortunate that the star formation 
rate density history has been mostly resolved to at least 
z < 6 , wi th the compilation data p resented in iHopkinsI 
(|2004l) and lHopkins fe Beacoml (|2006l) (see Figure Q] of this 
paper). Although unfortunately, the volumetric SN la 
rate history is much farther from consensus. While 

et al. 2010) show the promise of the large-scale surveys currently 
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Fig . 1. — The semi-analytical p*(z) mod el is show n (solid line), representing the compilation of measurements (grey points) from Hopkins 
2004). For compar ison, the best-f i tting the Cole ct al. (2001) parametric function (dashed line) is shown. This corrects an inaccurate repre- 
scntation shown in Valiantc ct al. (2009). Recent constrains from the HST+WFC3/IR Early Release data (Bouwens ct al. 2010; Ocsch et al. 
2010) are not included in the model fit. 



it is expected that the SN la rate sho uld demonstrate 
some increase with lookback time (e.g., iJorgensen et al.l 
Il997h . the scattered and often disagreeing rate measure- 
ments to date (most shown in Figure [5]) has made this 
evolution difficult to resolve. Q Figure [5] shows rate 
measure ments from various authors in several redshift 
regim es (ICappellaro et al.lll999t lReisdl2000t lHardin et all 



2000 1: iPain et all l2002t IStrolgerl 1200a iTonry et al.l 12 003; 



Madgwick et al.l 120031: [Blanc et all l2004t iBarris fe Tonryl 
2006; Poznanski et al. 2007; iKuznetsova et all [2008; 



Dah len et alll2008l: iDildav et alll2010D . It is presently un- 



clear where the point-to-point variation in the rate mea- 
sures stem from, although the most likely culprits are sur- 
vey completenesses (i.e., were all discoverable SNe la iden- 
tified?) , and differences in applied corrections for undis- 
covered events, typically referred to as control-time cor- 
rections. 

These discrepancies in SN la rate measures pose the 
greatest limitation in determining ^(t) empirica lly. The 
HST-SN survey rates (jDahlen et all |2004 |2008) have at 
least given a self-consistent measure of the SN la rates 
over a wide redshift range of 0.2 < z < 1.8 (see Figured]), 
making a determination for $(r) from these rate mea- 
sures alone reasonable. Additionally, the star formation 
rate density compilation data in th e z > 3 range is largely 
from the same GOODS/UD F fields (|Giavalisco et al.lfe OOl 
iBouwens et aDl2007l [20091 making the GOODS/HST-SN 



dataset most ideal for empirically probing the intrinsic 
$(r) function of SNe la, relatively free of the point-to- 
point Ria(z) biases, and biases due to cosmic variance. 

3. THE MODEL TESTS 

While the goal of this investigation i s to d etermine the 
$(r) function from the iDahlen et all (|2008l) data, a di- 
rect comparison to just the binned rate measurements is 
unsuitable as it provides only one or two datapoints to as- 
sess the production of SNe la at z > 1.2. A more robust 
comparison can be made to the individual SN la events 
that went in to the volumetric rate calculations. In gen- 
eral, the volumetric SN la rate is determined by comparing 
the survey SN la yield to the product of volume and the 
effective period (or control time) of the survey. In a red- 
shift range (Az = z 2 — z\) centered at a given redshift (z), 
the rate is generally expressed by: 



Ru(z) 



(4) 



E£ Z1 t> c (zi)]Av( z y 

where the Ni a (z) is the yield in the z— bin, t' c (z) is the 
control time, corrected for time dilation, and AV(z) is 
the conically sliced volume from z 2 to z\ (the bin size). 
A complete description of these parameters and how they 
are determi ned for the HST-SN s urvey, including t' c (z), is 
presented in lStrolger et all (|2004l ). 

We probe for an empirical $(t) function in a metho d 
also described more completely in IStrolger et all (|2004l ). 
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Fig. 2. — Various type la supernova rate measures as a function of redshift. Measures are averages over small rcdshift bins centered at the 
measured point. The widths of these bins are not shown for clarity in the diagram. 



To summarize, with assumed p+(z) and $(t) models, we 
use Equations [3] and H] to predict the expected redshift dis- 
tribution of SNe la for the survey. This result is compared 
to the observed redshifts of each SN la to produce a condi- 
tional probability test in an application of Bayes' method, 
where: 

P[p*(z), $(r)|Data] « P[Data|p*(z), $(r)], 

56 

=n N ia(^), 

1=1 

56 

= l[Ri a {z i )t' c (z i )AV(z i ). (5) 

The predicted number distribution, given the assumptions 
on the models, then serves as a probability function for 
finding SNe la at the specific redshifts in which we have 
found them. By this method, we maximize our leverage 
on the best $(t) model by using all 56 SNe la, rather than 
just the four binned rate measurements. 

We normalize the probability distributions to serve 
as a relative likelihood statistic. Changes in the input 
model parameters will allow changes in the likelihood with 
redshift. For this investigation, we use the extinction- 
corrected p*{z) determined from a semi-a nalytical fit to 
vario us measures (functional form shown in lStrolger et all 
l2004f) . The correction to extinction internal to each galaxy 
is consistently applied to both the published star forma- 
tion rate densities and to the control t i mes fo r the SN la 
rate calculations, both in iDahlen et al.l ()2008l ) and in this 



analysis. It is assumed that the p*(z) model and all other 
dependencies (e.g., Qm, ^a, Hq, and survey parameters) 
are sufficiently well determined that their uncertainties do 
not significantly contribute to the overall probability. 

The test method is to select a model $(r), calculate 
the Ni a (z), and determine a relative Bayesian likelihood 
of the chosen model from a comparison to the observed 
redshifts of our discovered SNe la. The likelihood values 
maximize when the predicted redshift distribution matches 
the observed redshifts, and minimize (or become zero) if 
for example a discovered SN la's redshift falls outside the 
predicted distribution. This high-risk test gives significant 
leverage in selecting the most likely $(r) when dealing 
with a statistically limited sample, free from biases associ- 
ated with binning data for \ 2 or K-S tests on observed red- 
shift distributions or rate measurements. Unfortunately, 
it also has a limited ability to address which $(r) models 
are implausible given the data, and therefore such inter- 
pret atk2ns_shcjfidbe made cautiously. 

In IStrolger et all (pOOl we examined a few simple 
single-parameter models (Gaussian models with variance 
tied to the mean, and exponential models), linearly iterat- 
ing through the characteristic delay time, f (describing the 
Gaussian mean time or the e-fold time), for a maximum 
likelihood. In that preliminary Bayesian test, the HST-SN 
data provided the highest likelihoods for Gaussian $(r) 
models with a mean of f « 3.5 Gyr and a T = 0.2, seem- 
ingly indicating a single, highly delayed mechanism for all 
SN la production. Models with incubation times less than 
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2.0 Gyr where inconsistent with the data to 95% confi- 
dence. However, the preliminary investigation was limited 
in that the tests could not assess large asymmetric skews 
(other than with an exponential decay), or the possibil- 
ity multi-modal populations. Most of the power of the 
previous tests was in determining the likeliest mean of in- 
cubation times, and there was limited ability to reveal a 
more complex $(t) distribution. 

As the natural next step in this investigation, we now 
test a more robust delay-time model, capable of more ac- 
curately reproducing the theoretical distributions for SD 
and DD models at one extreme, and (^-function delay times 
at the other. The unimodal, skew-normal $(t) function is 
defined as: 



*(r) 



exp 



-(r-Q ; 
2w 2 



exp 



dt', 



where location (£)□ scale (w 2 ), and shape (a) define the 
mode time (f, as defined in the previous tests), variance 
(cr 2 ), skewness (71), and kurtosis (72) of the model func- 
tion by, 



f = £ + lu5 



vTT^ 2 



= a/ 1 - 



7i 



28 2 



1(4-,) 

2 [ ; (l-2J 2 /7r) 3 / 2 ' 



72 = 2(tt - 3) 



(l-26 2 /7r) 5 



An illustration of the diversity of testable delay time dis- 
tribution functions from Equation [5] can be seen in Fig- 
ure |3l where the four 3>(t) models shown are created from 
independent choices of £, u>, and a model parameters. 

To test which model parameters values best fit our data, 
£, uj, and a have been jointly explored i n a Metropolis- 
Hastings Markov cha in Monte Carlo (jHastingsl Il970t 
Metropolis et al.lll953l hereafter MCMC) to find the most 
likely regions of the 3-D parameter space. To illustrate the 
MCMC test — a three-dimensional array was constructed 
for all £ (in range — 10 . . . + 10), v (range . . . + 10), and 
a (-10 . . . + 10), in intervals of 0.1. The MCMC values in 
this array were all initially set to zero. An initial starting 
point in £, u>, and a was randomly selected, and a relative 
Bayesian likelihood value was determined from Equation [5] 
at this starting point. The algorithm then randomly de- 
termined a step in all parameters simultaneously, where 
each step size (and direction) was determined from a nor- 
mal distribution, centered at zero with a — 0.25. In this 
way, 68% of the steps were A(£, lu, a) < ±0.25 (rounded to 
the 0.1 interval in the grid). The conditional test of Equa- 
tion[S]was then run on the new £ + A£, w+Aw, and a + Aa 
position. If the Bayesian likelihood value was greater at 
the new position than it was for the starting position (or 
equal to it), the algorithm incremented the MCMC array 



value in the new position by 1, adopted the new position 
as the starting point, and repeated the process of selecting 
a new step for evaluation. If, however, the Bayesian like- 
lihood for the new position was less than the value at the 
starting position, then the MCMC at the starting position 
was instead incremented by the ratio of likelihood values^ 
The starting position remained unchanged, and a new step 
was randomly selected and evaluated. 

After a "burn in" of 125 iterations (which were dis- 
carded), over 1000 test iterations were made where Equa- 
tion [5] served as the conditional likelihood between each 
step in the process. The MCMC array values thus built 
up at points in the grid with highest likelihood. Fig- 
ure [4] shows the most likely region of the MCMC test, 
in which > 95% of MCMC values lie. The peak of the 
confidence region is located at £ = 3.2+Q4, u = 0.2+gJ, 
and a — 2.2 ± 1.4 where the uncertainties are an approx- 
imation of the 95% confidence region. These model pa- 
rameters represent a $(r) with mode f = 3.4 Gyr and 
a = 0.14, shown in the upper right of Figure [4] The skew 
(71 = 0.51), and kurtosis (72 = 0.35), are nearly irrelevant 
for such a narrow distribution. The test indicates a $(t) 
that is surprisingl y consistent with th e best-fit narrow- 
gaussian model of Strolger et al.l ([2004) in that the mean 
incubation time is approximately 3.5 Gyr, and the entire 
distribution in contained within a 3 to 4 Gyr span. The 
implication is, once again, that there is only a single mech- 
anism for SN la production that requires almost exactly 
3.5 Gyr to develop WD progenitors into SNe la. The lack 
of width in the $(t) or a strong asymmetry is incorrob- 
orative of multiple components in this high redshift data. 
The SN la rate history, calculated from the MCMC best-fit 
delay-time distribution, is shown in Figure [5] The compar- 
ison rate models shown in the figure are scaled to match 
the observed SN la rates at z < 0.1. 

4. DISCUSSION 

The strong preference for a very specific delay time 
(with little variation) is uncanny, given the wide free- 
dom of the tested delay-time distri bution model. If 4>(t) 
we re intrinsically bimodal (as in M annucci et al. 2006 
or IScannapieco fc Bildster] I2005L for example), the test 
should have preferred a substantially wider and greatly 
skewed $(£,w,a) model, similar to what is shown in Fig- 
ure [3] A similarly exponential $(t) would have been ex- 
pected if the SN la rate were directly tied to just the avail- 
ability of WD at any given epoch, as is suggested by the 
log[$(r)] w —0.5 log(r) + constant model derived from 
a calculation of WD formation rates, and dominated by 
the m ain-sequence lifetime s of their 3 — 8 Mq progenitor 
stars (|Pritchet et al.ll2OO80 . But again, the MCMC test 
here does not find this to be the preferred model. The im- 
plication is that SN la are predominately characterized by 
a single channel (or mechanism) for explosion that requires 
almost exactly the same incubation time from formation 
to explosion. 

This interpretation, however, is difficult to reconcile 
with various results at lower redshifts. Studies in z < 
0.1 galaxies have shown strong evidence for SN la het- 
erogeneity that extends beyond the luminosity disper- 



2 Different from the initial mass function, §(M). 

3 The ratio is the Bayesian likelihood value (Equation [SJ at the new position over the likelihood at the starting position. 
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Fig. 3. — An illustration of the flexibility of the skew- normal model used in this analysis (from Equation |6)l . The g rey solid lines exem- 
plify $(r) models for diff erent choices of £, oj, and a, while the red lines are example model i 1 > (t) from the literature (Grcggio ct al. 2008, 
solid; Pritc het et al.H2008l . dotted; and lScannapieco fc Bildsten 2005, dashed). 



sion char acterized by t he luminosity-hghtc urve width 
relations (|Phillipsl Il99l iPhillips et all 1 1999( 1. including 
low-luminosi ty SNe la with strong Ti absorption (e.g., 
SN 1991bg. [Ruiz-Lapuente et all[l9 93;l high- luminosity 
SNe la wit h weak Si II absorptio n at early times (e.g., 
SN 1991T. iFilippenko et al.lll992l) . and several examples 
of individually peculiar supernovae. In addition, rate mea- 
sures in low-z galaxies show that late- type galaxies are 
rather prolific producers of SNe la, producing an order 
of magnitude more events than early-type galax ies of the 
same total stellar mass (Ma nnucci et al.l [2005). More- 
over, trends show that the most luminous SNe la are 
largely absent in early-type galaxies, and l ow- luminosity 
event s are deficient in lat e-type galaxies (|Hamuv et al.l 
I2000t lAltavilla et~aTI I2004D . The implication has been 
that the rate of recent star-formation has a significant im- 
pact on the production of SNe la, and that the overall 
rate o f SNe la is a sum of prompt and delayed compo- 
nents. IScannapieco fe Bildsten (2005)) have argued that a 
two-component model with a large contribution of prompt 
SNe la (20-40% of all SNe la) can resolve the inter-cluster 
Fe content and [O/Fe] abundance in the Milky Way and re- 
main consistent with the observed low-z rates. Moreover, 
the rate at higher redshifts (z ~ 0.6) from the Supernova 
Legacy Survey also support the two-component model, but 
with a prompt component increased in its f raction to as 
high as 80 - 90% by z < 1 dNeill et alll2006h . 



Indeed, a more direct relationship between SN la rates 
and star formation rates (with little-to-no delay) appears 
possible given notable similarity between the p*(i) func- 
tion and most Ri a (i) measurements at z < 1 in Figure (5[ 
using a scaling where approximately one SN la is produced 
for every 500 M Q created with essentially zero delayH It 
is, therefore, interesting to address what results when the 
HST-SN sample is cut to z < 1.0 and our MCMC test are 
re-done. 

4.1. Cuts on the Data: The z < 1 Sample 

We performed our MCMC test on a sub-sample of 
the HST-SN survey SNe la, selecting only the 32 events 
with z < 1.0. The resulting MCMC likelihood con- 
tours are shown in Figures [6] and [7J The results show 
the highest likelihoods for power-law <1>(t) models which 
show long tails to large delay times (peak at £ = 
0.6, w = 7.5, a = —1.0). These models bare some 
resembla nce in shape to several S D and DP models, 
includ i ng lYungelson fe Liviol ([2000), M atteucci fe Recchl 
(I2001h. iGreggiol (120051) . " and iGreggio et al.l (|2008h . and 
the iPritchet et all (|2008l) WD availability model. The 
match is not perfect as literature models are typically 
steeper, requiring disproportionately more prompt (< 
lGyr) SNe la than the results of our MCMC test. The mis- 
match may be partly due to some rigidity in skew-normal 
model of our MCMC test, in that the $(£, lu, a) model may 



This is effectively reduces Equation \3\ to Ri a (t) = ep*(t), where e = 0.002 M, 
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Fig. 4. — Results of the Markov chain Monte Carlo test of three-parameter skew-normal distribution, where £, a), and a (which describe 
the shape of the delay time function) are explored. The best-fit <J>(r) has be en determined from centroid of the 95% likelihood region, and 
is shown in in the upper-right panel (s o lid lin e). Shown for comparison is the Grcggio ct al. (2008) single degenerate <E>(t) model (G08; solid 
grey line), the [Scannapicco & Bildstcn (2005) bimodal model (SB05; dashed grey line), and the lPritchet et all (120081 ) white-dwarf availability 
model (P08; dotted grey line). The scalings of the < J > (r) models are arbitrary. 



not be flexible enough to reproduce broken power-law-like 
distributions. However, tests of the flexibility of the model 
(exemplified in Figure s uggest this is not the case. 

A recent investigation by I Raskin et al.l (|2009al) compar- 
ing the locations of SN la with the distribution of light in 
z < 0.07 spiral ga laxies (testing the p aren t population of 
stars , as is done in lFruchter et aLlfeoOd and lSvensson et al.l 
20f for Gamma-ray Bursts) offers an alternative explana- 
tion. Their results seem to indicate a substantial delay for 
even the potentially prompt SN la population of at least 
200-500 Myr. Although this is not as large as the delay 
from the MCMC results shown in ^ a scenario in which 
events delayed by < 0.5 — 1 Gyr simply do not occur (es- 
sentially cutting all delays below 1 Gyr in the upper-right 
panel of Figure [6]) would improve the agreement of our 
z < 1 sample results with many $(r) models. 

ft is also interesting to note that the z < 1.0 MCMC 
results show a second slightly smaller likelihood peak near 
f = 4.2 Gyr (£ = 4.5, cj = 0.5, a = -1.0; see Fig- 
ure [7|) . It is tempting to interpret the existence of two 
likelihood peaks as support for bimodality in the intrin- 
sic delay-time distribution, implying two separate SN la 
mechanisms. However the more cautious interpretation is 
that this MCMC test is not designed to find bimodal $(r) 
models. Its design is to find the peak and distribution for a 
model that has a single mode, assuming the intrinsic $(r) 
can be accurately characterized in this way. Secondly, the 



result also only materializes when a sub-sample of z < 1 
data is used, and not with the full dataset. There is cur- 
rently no reason to suspect that our SNe la at z < 1 carry 
more weight than those at z > 1, although the possibility 
is discussed further in £14.31 The MCMC result from the 
full sample should be principally considered. 

4.2. Cut on SN la Peak Luminosity 

Another intriguing cut on our sample is to see if the 
redshift distribution is dependent on event luminosity. 
The Supernova Legacy Survey noted an interesting trend 
in which their lowest redshift sample (0.0 < z < 0.1) 
were dominated by lower luminosity (narrower lightcurve 
width) events, and their highest red shift sample (0-75 < 
z < 1.5) by more luminous events (jHowell et al.l l2007h . 
If low-luminsoity events have a substantial fraction of 
SN 1991bg-like events, and these events represent a sepa- 
rate channel for SN la production (perhaps more delayed 
than for normal SNe la), then this subset could show a 
different delay-time distribution than the parent popula- 
tion, reflected in the observed distribution of redshifts for 
these events. The converse could also be expected for the 
high-luminosity sample, which should contain a high frac- 
tion of SN 1991T-like events that could also represent a 
separate SN la mechanism. A caveat is that by nature 
less luminous events will be less numerous in the highest 
redshift regimes, as magnitude-limited surveys lose sensi- 
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Fig. 5. — Model type la supernova rates compared to rate measurements, displayed in terms of the age of the universe (concordance 
model cosmology is assumed). Upper panel: The solid black line is the best-fitting R/ a (t) from the 3>(t) model deri ved in this inv estiga- 
tion, with an effective mea n r = 3.5 Gy r. For comparison, the blue li nes show model Hj a (t) using 3>(t) models from Grcggio ct al. (2008, 
solid), Pritchet ct al. (2008, dotted), and Scannapicco & Bildstcn (2005, dashed). The red line (and secondary axis) is the model p*(t). Lower 
panel: The best-fit model shown with a confidence region (dotted black lines). 



mode near 3.2 Gyr) than the less luminous sample (with a 
mode near 4.2 Gyr and a broader distribution). The lack of 
great change is expected, due to appropriate corrections in 
the control time calculations which truncate the assumed 
luminosity functional at high or low luminosity. 

4.3. Sensitivity of the HST-SN Survey 

Arguably, most of the rationale for such large delay 
times comes from the notable dearth of SNe la z > 1.4 
from the survey. There has been some concern on the 
completeness and sensitivity of the HST-SN survey to 
events in the highest magnitude (and redshift) ranges, 
but recent attempts find more SNe la at z > 1.4 in 
the same HST-SN data via independent detection crite- 
ria (|Kuznetsova et all 2008f). novel pixel- bv-pixel N(N — 
l)/2 comparisons ( Rodney fc To nrv 2007), and a thorough 
revie w of the deepest region of the GOODS South (the 
UDF; IStrolger fc Riessl [2006D have failed to produce con- 
vincing additional candidates. Additionally, there is fur- 
ther evidence for a lack of z > 1.4 SN e la from the Subaru 
Deep Survey (jPoznanski et al.ll2007T) . 

Ano ther critical concern discussed in iGreggio et al.l 
( 2008) is that high redshift galaxies are likely to have more 
internal extinction than their low-z counterparts due to 

5 AMi 5 (B) = f .07 is the mode of SN la lightcurve widths {Phillips et al.llT999l '). 

° A gaussian distribution consistent with recent results from the Sloan Digital Sky Survey l lYasuda fc Fuku gita 2010). 



tivity to them at slightly lower- z (our survey sensitivity is 
discussed further in £)4.3[) . 

Using fits to the lightcurve widths, determined in the 
rest-frame of each SN, we separate the sample into high- 
luminsoity (with AMi5(B) < 1.07) and low-luminosity 
(AM 15 (B) > 1.07) events^ It should be noted that with 
the exception of the 22 events used for the cosmological 
investigation (|Riess et al.ll2007D . most of the remaining su- 
pcrnovae have limited lightcurve information, which limits 
further precise determination of peak luminosity or decline 
rates. Figure [5] shows the redshift distribution for the high 
and low luminosity samples, compared to the full sarn- 
ie. A s expected, there is a trend similar to iHowell et al.l 
2007h for the less luminous supernovae to shift to lower 
redshifts, peaking near z ~ 0.8 rather than z ~ 1.0, and an 
opposite trend (although less pronounced) in the more lu- 
minous sample. These shifts, however, do not translate to 
significant changes in the derived delay-time distribution 
function. 

We performed our MCMC test independently on the 
bright (AMi 5 (B) < 1.07) sample of 26 events, and the 
faint (AMi5(B) > 1.07) sample of 30 events, the results of 
which are shown in Figure [5] As is shown, the more lumi- 
nous events do prefer slightly shorter delay times (with a 
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Fig. 6. — Same as for Figure[4]but on the sub-sample of HST SNe la with z < 1.0. 




the very enhanced dust production associated with the 
high rate of star formation. At these redshifts, we probe 
the rest-frame UV region of the SN la spectrum, a region 
which is inherently photon-deficient, and extremely sen- 
sitive to host extinction. Naturally, attempts are made 
to correct for the portion of SNe la lost to internal ex- 
tinction through the control times (essentially "efficiency 
corrections"), by using modeled and observed radial dis- 
tributions of SNe la i n low-z galaxies and extinction distri - 
butions within them (jHatano et al.lll99&|j"ha et alJll999l) . 
However this concern is for an additional extinction over 
that which is traditionally account ed for. To date, nei- 
ther the HST-SN cosmology data (|Riess et al.l [20071 nor 
programs which span intermediate redshift ranges (e.g., 
SNLS and ESSENCE) and bridge the gap from low-z to 
high-z, show any evidence for an "extinction excess" trend 
with redshift. There is, in fact, an opposite tendency for 
routines such as MLCS2k2 (which fit Rv as a free pa- 
rameter in lightcurvc fitting) to prefer "grayer" extinction 
laws for higher redshift SNe la, with 1.6 < (Rv) < 2.7 
pha. Riess. fe Kirshnerll2007ft . 

4.4. Changes in the Star- Formation Rate Density, 
or a Possible Metallicity Effect on SN la Production 

Our best-fit ^(t) model is very different than what may 
have been expected from a more classical modeling of 
the SN la progenitors, and from recent interpretations of 
lower- z data. It is difficult to see how both could be cor- 
rect. However, there is a way out of this apparent quag- 
mire. If indeed the standard SD and DD $(t) models are 



correct, then there must be an additional factor which puts 
additional constraint on the shape of the derived $(r) at 
low t in this study. 

One very likely possibility is that the shape of the p+(z) 
function (shown in Figure[T]) does not flatten beyond z > 6, 
rather declines more sharply. Very recent (and tentative) 
results from the HST+WFC3/IR early release portion of 
an 192-orbit ultra-deep survey of the HUDF (Illingworth, 
GO 11 563) seem to indi cate that p+(z) declines to nearly 
z = values near z = 8 dBouwens et al.l 2010t lOesch et al.l 
l2010fh Based on test done in lStrolger et al.l (|2004l) on al- 
ternative models for the p*(z) function (including a model 
with a sharper decline), we predict that similar MCMC 
test will show greater likelihood values in smaller r re- 
gions, but the best-fit overall will remain high. 

Another potential culprit could be an innate metallic- 
ity effect, which suppresses the prompt component of the 
natural $(t) through some requirement of a "minimum 
metallically" for the SN la mechanism. A physical ratio- 
nale for a minimum metallicity effect in SD scenarios could 
be that potential WDs progenitors must develop a suffi- 
cient counter-wind to allow for steady mass accretion, to 
bulk up the core of the WD without triggering surface H 
& He flashes (novae) and substantial mass losses, or accre- 
tion induced core-collapse supernovae, o r other scenarios 
which ultimately fail to produce a SN la ([Kobavashi et al.l 
1 1998k iKobavashi fc Nomotoll2007t ) . This is s omewhat sup- 
portiv e of the resent observational results of ICooper et al.l 
( 2009]) on the large-scale environmental impacts on SN la 
production, but acts in the opposite way than what would 
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Fig. 7. — Same as for Figure|6] The axes have been expanded to reveal the second likelihood peak. 



be expected from their results. 

If metallicity is a critical missing factor, a new parame- 
terization of the^SN la rate with redshift could then be: 

R/«(t)=e / $(t-t')p*(t')C{[0/H](t-t')}<ii', (7) 

Jt 

Where C{[0/H](r)} describes the efficiency in success- 
fully making SNe based on the [O/H] in the system at 
the time of formation. The [O/H] ratio is chosen as it 
is suspected that products of the CNO-cycle, specifically 
22 Ne which is traced by 16 abundances in the interstel- 
lar medium, have a larger impact o n the SN la outcom e 
than Fe-peak elements in the ISM (jTimmes et al.ll2003l) . 
although it is not at all clear how ISM metallicity impacts 
SN la luminosity or production. The global metallicity 
enrichment of the universe sh ould be a slowly decre asing 
function with lookback time (jKulkarni et al.l |2005[ ) , but 
due to its convolution with the intrinsic delay-time func- 
tion of SNe la, could pr ovide a steep cutoff to the SN la 
production at highest- z (|Riess fc Liviol l2006'). 

5. CONCLUSIONS 

The HST-SN data are most supportive of a single dom- 
inant mechanism for the production of SN la, which re- 
quires between 3 and 4 Gyr of incubation from system 
formation to explosion. This mechanism is supportive of 
single degenerate models in which the companion donor 
stars are low mass (< 2 M Q ), based on the main-sequence 
lifetimes which dominate (much more so than the accre- 
tion times) the incubation period. The data are largely in- 
consistent with progenitor scenarios with short (< 1 Gyr) 



development times, and are inconsistent global scenarios 
where prompt mechanisms make up a substantial fraction 
of all channels employed by SN la progenitors. 

The preference of our delay time model tests for high t 
are largely motivated by the observed reduction in z > 1 
SNe la. Further investigations of the SN la rate at even 
higher redshifts, from 1.5 < z < 3.0, should elucidate 
possible metallicity trends (or other effects) from sensi- 
tivity issues, and are plausible in large campaigns with 
HST with the IR channel of WFC3. They will also be 
easily assessable in future space-based missions such as 
the James Webb Space Telescope. In the meanwhile, it 
is also imperative that the SN la rate measures in the 
0.1 < z < 1.0 range reach consensus to make more rigor- 
ous global comparisons to the SN la rate history, and fur- 
ther refine the empirical delay-time distribution function. 
This is within reach with large-scale surveys such as the 
Palomar Transient Factory and Pan-STARRS. Ultimate 
comparisons will be achievable when the Large Synoptic 
Survey Telescope, and the NASA/DOE Joint Dark Energy 
Mission/International Dark Energy Cosmology Survey are 
realized. 
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Fig. 8. — Redshift distribution of SNe la from the HST-SN survey (black histogram). The red histogram shows the distribution for 
AMi5(B) > 1.07, and the blue histogram is for AMis(B) < 1.07, in the rest-frame of each supernova. 
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AMis(i?) > 1.07 sub-sample (right four panels). 



and the less luminous 



cial support for this work was provided by the Kentucky lege and Fellowship Program, and by the Western Ken- 
Space Grant Consortium (K07R10), funded by a NASA tucky University Research Foundation. 
Training Grant as part of the National Space Grant Col- 



REFERENCES 



Altavilla, G., et al. 2004, MNRAS, 349, 1344 
Barris, B. J., & Tonry, J. L. 2006, ApJ, 637, 427 
Blanc, G., et al. 2004, A&A, 423, 881 

Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, ApJ, 
670, 928 

Bouwens, R. J., et al. 2009, ApJ, 705, 936 



Bouwens, R. J., et al. 2010, ApJ, 709, L133 

Cappcllaro, E.. Evans, R., & Turatto, M. 1999, A&A, 351, 459 

Cole, S., et al. 2001, MNRAS, 326, 255 

Cooper, M. C, Newman, J. A., & Yan, R. 2009, ApJ, 704, 68 
Dahlcn, T., et al. 2004, ApJ, 613, 189 

Dahlcn. T., Strolger, L.-G., & Riess, A. G. 2008, ApJ, 681, 462 



12 



Strolger , Dahlen, & Riess 



Dave, R. 2008, MNRAS, 385, 147 
Dilday, B., et al. 2010, arXiv:1001.499 
Filippenko, A. V., et al. 1992, ApJ, 384, L15 
Fruchter, A. S., et al. 2006, Nature, 441, 463 

Gallagher, J. S., Garnavich, P. M., Caldwell, N., Kirshner, R. P., 
Jha, S. W., Li, W., Ganeshalingam, M., & Filippenko, A. V. 2008, 
ApJ, 685, 752 

Giavalisco, M., et al. 2004, ApJ, 600, L93 

Grcggio, L. 2005, A&A, 441, 1055 

Greggio, L., Renzini, A., & Daddi, E. 2008, MNRAS, 388, 829 
Hamuy, M., Trager, S. C, Pinto, P. A., Phillips, M. M., Schommer, 

R. A., Ivanov, V., & Suntzeff, N. B. 2000, AJ, 120, 1479 
Han, Z., & Podsiadlowski, P. 2008, IAU Symposium, 252, 349 
Hardin, D., et al. 2000, A&A, 362, 419 
Hastings, W. K. 1970. Biometrika, 57, 97 
Hatano, K., Branch, D., & Deaton, J. 1998, ApJ, 502, 177 
Hopkins, A. M. 2004, ApJ, 615, 209 
Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 
Howell, D. A., Sullivan, M., Conley, A., & Carlberg, R. 2007, ApJ, 

667, L37 

Howell, D. A., et al. 2009, ApJ, 691, 661 
Jha, S., et al. 1999, ApJS, 125, 73 

Jha, S., Riess, A. G., & Kirshner, R. P. 2007, ApJ, 659, 122 
Jorgensen, H. E., Lipunov, V. M., Panchenko, I. E., Postnov, K. A., 

& Prokhorov, M. E. 1997, ApJ, 486, 110 
Kobayashi, C, Tsujimoto, T., Nomoto, K., Hachisu, I., k Kato, M. 

1998, ApJ, 503, L155 

Kobayashi, C, & Nomoto, K. 2007, arXiv:0801.0215 

Kulkarni, V. P., Fall, S. M., Lauroesch, J. T., York, D. G., Welty, 

D. E., Khare, P., & Truran, J. W. 2005, ApJ, 618, 68 
Kuznetsova, N., et al. 2008, ApJ, 673, 981 

Madgwick, D. S., Hewett, P. C, Mortlock, D. J., k Wang, L. 2003. 
ApJ, 599, L33 

Mannucci, F., Delia Valle, M., Panagia, N., Cappellaro, E., Cresci, 
G, Maiolino, R., Petrosian, A., k Turatto, M. 2005, A&A, 433, 
807 

Mannucci, F., Delia Valle, M., k Panagia, N. 2006, MNRAS, 370, 
773 



Matteucci, F., k Recchi, S. 2001, ApJ, 558, 351 

Metropolis, N., Rosenbluth, A., Rosenbluth, M., Teller. A., and 

Teller, E. 1953, Journal of Chemical Physics, 21, 1087 
Neill, J. D., et al. 2006, AJ, 132, 1126 
Oesch, P. A., et al. 2010, ApJ, 709, L16 
Pain, R., et al. 2002, ApJ, 577, 120 
Phillips, M. M. 1993, ApJ, 413, L105 

Phillips, M. M., Lira, P., Suntzeff, N. B., Schommer, R. A., Hamuy, 

M., k Maza, J. 1999, AJ, 118, 1766 
Poznanski, D., et al. 2007, MNRAS, 382, 1169 

Pritchet, C. J., Howell, D. A., k Sullivan, M. 2008, ApJ, 683, L25 
Raskin, C, Scannapieco, E., Rhoads, J., k Delia Valle, M. 2009. 
ApJ, 707, 74 

Raskin, C, Timmes, F. X., Scannapieco, E., Diehl, S., k Fryer, C. 

2009, MNRAS, 399, L156 
Reiss, D. 2000, Ph.D. Thesis, University of Washington 
Riess, A. G, et al. 2004, ApJ, 600, L163 
Riess, A. G, & Livio, M. 2006, ApJ, 648, 884 
Riess, A. G, et al. 2007, ApJ, 659, 98 

Rodney, S. A., k Tonry, J. L. 2007, Bulletin of the American 

Astronomical Society, 38, 807 
Rosswog, S., Kasen, D., Guillochon, J., k Ramirez-Ruiz, E. 2009. 

ApJ, 705, L128 
Ruiz-Lapuente, P., et al. 1993, Nature, 365, 728 
Scannapieco, E., k Bildsten. L. 2005, ApJ, 629, L85 
Strolger, L.-G. 2003, Ph.D. Thesis, University of Michigan 
Strolger, L.-G., et al. 2004, ApJ, 613, 200 
Strolger, L.-G., k Riess, A. G. 2006, AJ, 131, 1629 
Svensson, K. M., Lev an, A. J., Tanvir , N. R., Fruchter, A. S., k 

Strolger, L. -. 2010. arXiv:1001.5042 
Timmes, F. X., Brown, E. F., k Truran, J. W. 2003, ApJ, 590, L83 
Tonry, J. L., et al. 2003, ApJ, 594, 1 

Valiante, R., Matteucci, F., Recchi, S., k Calura, F. 2009, New 

Astronomy, 14, 638 
van Dokkum, P. G. 2008, ApJ, 674, 29 
Yasuda, N., k Fukugita, M. 2010, AJ, 139, 39 
Yungelson, L. R., & Livio, M. 2000, ApJ, 528, 108 



